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We present a study of the elastic exciton-electron (X — e~) and exciton-hole (X — h) scattering 
processes in semiconductor quantum wells, including fermion exchange effects. The balance between 
the exciton and the free carrier populations within the electron-hole plasma is discussed in terms 
of ionization degree in the nondegenerate regime. Assuming a two-dimensional Coulomb potential 
statically screened by the free carrier gas, we apply the variable phase method to obtain the excitonic 
wavefunctions, which we use to calculate the Is exciton-free carrier matrix elements that describe 
the scattering of excitons into the light cone where they can radiatively recombine. The photon 
emission rates due to the carrier-assisted exciton recombination in semiconductor quantum-wells 
(QWs) at room temperature and in a low density regime are obtained from Fermi's golden rule, 
and studied for mid-gap and wide-gap materials. The quantitative comparison of the direct and 
exchange terms of the scattering matrix elements shows that fermion exchange is the dominant 
mechanism of the exciton-carrier scattering process. This is confirmed by our analysis of the rates 
of photon emission induced by electron-assisted and hole-assisted exciton recombinations. 

PACS numbers: 71.35.-y, 78.55.-m, 78.55.Cr, 78.55.Et, 78.67.De 

I. INTRODUCTION 

Scattering processes involving excitons are at the heart of a wide range of phenomena in semiconductor optics. 
Excitons are often treated as elementary bosons, but the complexity of the scattering problem obviously lies in their 
composite nature: excitons are Coulomb-correlated quasiparticles made of two fermions, a conduction band electron 
and a valence band hole. It is only recently that a many-body theory of exciton scattering was put on firm grounds in 
a series of papers by Combescot et al (see e.g. Refs. [H, and references therein). Their formalism, based on fermion 
commutation techniques, allows one to calculate correctly the scattering matrix elements as well as the transition rates 
of two excitons. In the present paper, we are only concerned with the simpler problem of exciton-carrier scattering, i.e. 
a 3-body problem with well defined interaction potentials between the three scattering partners, including exchange 
effects explicitly. More precisely, we are interested in exciton-free carrier scattering in III-V GaAs-based and II- VI 
ZnSe-based quantum-well systems. 

Early investigations of the microscopic mechanisms yielding low lasing thresholds and high optical gain in bulk 
semiconductor lasers included detailed theoretical and experimental studies of radiative recombination involving ex- 
citon scattering. Various scattering processes including X — X, X— LO-phonon, X — e~ scattering [!, 0, HI and 
X — h scattering 0, [|| were studied and compared, and it was reported that X — e~ scattering is a process yielding 
important optical gain in bulk materials [H, |j,[l|. In particular, Haug and co-workers gave a detailed study of the 
temperature dependence of the lasing thresholds for each of the scattering process [4|, [5j. Later, Feng and Spector 
studied the exciton-free carrier elastic and inelastic scatterings in quantum wells to compute the related cross sections 
and exciton linewidths, but they did not account for carrier exchange in their analysis ,6|. Here, the aim of our work is 
two-fold: first we want to compare the respective contributions of the direct and the exchange terms to the scattering 
matrix elements; then we want to study quantitatively the exchange term dependence on the carrier effective mass 
and its impact on carrier-assisted exciton radiative recombination in quantum wells. 

The interest for extensive studies of photoluminescence (PL) spectra @, 1, i, El, El El, El El, El, El, E3, El, 
El, HI, [U, H3 lies in the fact that they are useful to study a rich variety of phenomena in semiconductor physics 
and allow non-destructive characterization of semiconductors. In the case of mixed exciton/electron-hole plasmas, 
many processes occur that contribute to PL spectra. These processes include: direct exciton recombination, electron- 
hole recombination and several scattering processes involving quasi-particles such as excitons, phonons and excited 
carriers. The importance of each contribution to PL depends on the thermodynamic properties of the plasma and 
hence on the excitation conditions as well as the type of material under consideration that determine the populations 
and temperatures of each quasi-particle gas. Moreover, the dimensionality of the system is also of great importance 
since quantum confinement not only enhances the effective strength of the Coulomb interaction between the carriers, 
but changes drastically the density of states of the particles. 
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Early works by Kubo, Martin and Schwinger (KMS) [23j, [24J suggest that in quasi-equilibrium the intensity of PL is 
proportional to the absorption coefficient times the Bose distribution describing the photon gas interacting with the 
semiconductor medium. Detailed studies of the excitonic resonance features and their dynamics through PL spectra 
based on the KMS approach led to the interpretation of the build-up of the excitonic resonance below the gap as direct 
evidence of excitonic formation on the sub-nanosecond timescale, see e.g. Refs. SHI 13. However, calculations 
by Kira et al [l?], EH clearly show that there is no straightforward connection between exciton formation and the 
presence of an excitonic peak in the emission: free plasma emission can occur at the exciton energy. Galbraith et 
al [l9j reported later on the existence of luminescence at the exciton energy in GaAs-based multiple quantum well, 
consistent with the theoretical work of Kira et al [l7], Ell, and Chatterjee et al [2(| identified conditions under which 
the PL emission at the exciton resonance after excitation in the continuum is dominated by the free electron-hole 
plasma or by an incoherent exciton population at low temperature (but they could not determine experimentally the 
fraction of excitons contributing to PL). This contradicts the findings reported by Szczytko et al [2l| who claimed 
that excitons provide the dominant contribution to the luminescence signals at the exciton energy, and concluded 
that for densities, temperatures, and time scales actually used in time- resolved experiments the Coulomb correlated 
plasma contribution may be neglected. The conclusion on the role of a finite exciton population based on the results 
obtained from simple rate equations used by Szczytko et al may require further investigation to be confirmed, but 
models that take into account the existence of a finite population of excitons in excited QWs deserve nonetheless 
attention. 

In our model, we assume that the exciton and carrier populations at room temperature in a low density regime have 
reached thermodynamic equilibrium. During the course of the scattering process, an exciton leaves its initial state 
and reaches the photon line before radiative recombination, transferring its momentum to its scattering partner, a 
free carrier. The exciton-electron scattering was found to be an important mechanism producing photoluminescence 
in bulk materials i, i, S El and optical gain in thin films [22j • It is thus of interest to investigate this mechanism 
in mid-gap and wide-gap quantum wells, and also obtain further insight by including the exciton-hole scattering in 
the analysis. Therefore, the object of the present work is not the calculation of the full PL spectra for which a full 
microscopic theory [I?], 0, HcJ is required, but the quantitative study of the relative importance of X — e~ and X — h 
scattering mechanisms, including exchange. Direct exciton recombination, known to yield a strong emission at exciton 
energy QQ, [ll[ , as well as scattering processes other than that of excitons with free carriers are beyond the scope of 
the present work and hence are not included in the paper. 

A complete treatment of the balance between bound excitonic states, scattering states and free carriers is a difficult 
problem that has attracted attention for several years [ID, H|| [13, [H, [H, [3(| HU . In the case of QWs, advances were 
made in Ref. [3(| where the degree of ionization, a, of a non-degenerate two-dimensional (2D) electron-hole plasma 
was calculated as the ratio n°/n of the density of free carriers, n°, to the total plasma density, n — n° + n co " where 
n co " is the density of Coulomb correlated carriers [2(|. To account for occupation of scattering states and overlap 
of exciton wavefunctions at moderate and high densities this concept was further refined in Ref. [3l[, so that the 
definition of a is given by a = 1 — n a t>/"-a, where n a is the total density of carriers of type a and n a b originates 
from the interaction of carriers of type a with carriers of type b (a ^ b). The ionization degree is now formulated 
in terms of elementary excitations (electrons and holes). Note that in the low density limit, n a b can be identified 
to the exciton density. The deg ree of ionization whose values are between and 1, is both density and temperature 
dependent. Authors of Ref. '30J demonstrated that for wide-gap semiconductor quantum-wells at room temperature, 
the equilibrium consists of an almost equal mixture of correlated electron-hole pairs and uncorrelated free carriers 
whereas this is not the case for mid-gap materials. 

For a plasma temperature of 300 K, we ensure the occupation density is small enough (nondegenerate limit) to 
neglect phase-space filling effects. Under these conditions, the scattering between excitons and carriers involves mostly 
thermalized Is excitons with a finite center of mass momentum and free carriers. This assumption is valid at room 
temperature and moderate densities where the Is exciton population dominates all other Coulomb correlations in the 
electron-hole plasma. Calculations in an unscreened Coulomb potential with hydrogen atom-like exciton wavefunctions 
have been reported by Kavokin [32]. In the present work, the electron-hole bound states are computed assuming a 
statically screened Coulomb potential (whose asymptotic behavior permits the application of the variable phase 
method [Hj]). 

In Section II we detail the basic assumptions made for the model of the interacting 2D electron-hole plasma and 
define the ionization degree. Section III is devoted to the derivation of the exciton wavefunctions from the variable 
phase theory. In Section IV we calculate the exciton-free carrier scattering matrix elements that we include into 
Fermi's golden rule to compute the related emission rates. Numerical results are presented and discussed in Section 
V, where photon emission rates induced by exciton-carrier scattering in ZnSe and GaAs QWs are compared. 
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II. INGREDIENTS OF THE MODEL 



A. Underlying assumptions 

In this work we assume the quantum- well structures to be ideal two-dimensional systems in a two-band model with 
parabolic dispersion, neglecting detail related to valence-band mixing. We consider a neutral low density, 2D plasma 
composed of interacting electrons and holes, at room temperature. We thereby remain in the Boltzmann regime 
denned by 0: 



n\ 2 m Jg<^l, (1) 

where n is the 2D carrier density, A„ lc = (2ttH 2 / rUcksT) 1 / 2 the thermal wavelength and g the spin degeneracy of the 
confined carriers of effective mass m c . For T = 300 K, Eq. (TTJ) is satisfied for n <C 1.7 x 10 12 cm 2 for ZnSe-based QWs 
and n«r.2x 10 11 cm 2 for GaAs-based QWs [U- 

Due to spatial confinement of the carriers in QWs the Coulombic interaction between electrons and holes is enhanced. 
However, even at moderate plasma densities carrier-carrier interactions are weakened because of the density-dependent 
screening. Screening is one of the most important manifestations of the complex many-body interaction in the electron- 
hole plasma and the simplest approach is the use of the 2D statically screened potential [3^ : 



v s ( P ) = e - r^a q , (2) 

where Jo is the Bessel function; q s , in the Boltzmann limit, is known as the 2D Debye-Hiickel screening wavenumber 
(which depends on temperature and carrier density) , e is the static dielectric constant of the semiconductor and p the 
inter-carrier distance. Equation ([2]) describes the electron-electron and hole-hole repulsion; the attractive electron-hole 
potential is obtained by changing the overall sign. 

The above integral may also be written in terms of special functions: 



V s (p) = - 
e 



~-J (H (q s p) - N {q s p)) 
P 2 



(3) 



where Nq and Ho arc the Neumann and Struve functions respectively [36l , l37j . 

The total screening is the sum of electron and hole plasma screenings. As such we neglect the weak screening by 
neutral excitons. 



B. The degree of ionization 

In QWs where Coulomb forces are enhanced because of the spatial confinement of the electrons and holes, we 
expect a finite population of Coulomb correlated quasi-particles in the plasma even at room temperature. The 
balance between those correlations and the free carriers is obtained by calculating the degree of ionization, a [3~fl ] : 

n c h n e h /.x 

a = 1 = 1 . (4) 

n e nil 

In the above equation, the total density of a carrier of type a is defined as [31]: 



ri a = n° a + n aa + n ahl (5) 

where n a is the density of free carriers of type a and n aa originates from the interaction between carriers of the same 
type (ee or hh) . When a is close to unity, the thermodynamic properties of the electron- hole plasma are those of the 
ideal gas (defined by a — 1). For lower values of the degree of ionization the thermodynamic properties of the plasma 
deviate from those of the ideal gas and Coulomb correlations have to be considered. 

Detailed calculations of the degree of ionization and further discussion on the statistical mechanics of the 2D 
electron-hole plasma can be found in Refs. [29|, [3(J HH . Here we only stress that the degree of ionization is obtained 
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from the calculation of the partition function of the 2D electron-hole plasma. Scattering state contributions to 
the partition function have to be considered in addition to the bound-state sum as a proper account of scattering 
eliminates singularities in thermodynamic properties of the non-ideal 2D gas caused by the emergence of additional 
bound states as the strength of the attractive potential is increased. Inclusion of the scattering states also leads to 
a strong deviation from the standard law of mass action. Additionally, note that exchange between carriers within 
the electron-hole plasma is taken into account in the calculation of the electron-electron and hole-hole parts of the 
partition function [301 ]. 

The ionization degrees, acaAs and aznSe, are shown as functions of plasma density n on Fig. [1] for GaAs and 
ZnSe quantum- wells at temperature T = 300 K. In the nondegenerate regime, the function aGaAs(n) is monotonic 
whereas a-z n Se( n ) exhibits a minimum at low density. The presence of the minimum is due to the influence of Coulomb 
screening on the density dependence of the ionization degree: at very low density screening is negligible and n™ rr is 
proportional to the square of the total plasma density so a is a decreasing function of n; however with n increasing 
further, plasma screening becomes more and more important and hence reduces the strength of Coulomb correlations 
so nl°" is a decreasing function of n and a increases. This density dependence is more dramatic in ZnSe QWs where 
Coulomb forces are much stronger than they are in GaAs QWs. Note that aGaAs(^) and aznSe(«) not only differ 
qualitatively but also quantitatively: in GaAs QWs the electron-hole plasma is dominated by free carriers whereas it 
contains comparable densities of correlations and free carriers in ZnSe QWs. Knowledge of the degree of ionization is 
of great importance as the emission mechanisms depend on the nature of the plasma. For further detail see Ref. [30| . 
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FIG. 1: Ionization degree evaluated for GaAs and ZnSe quantum wells in the Boltzmann regime for T = 300 K, as a function 
of the 2D plasma density scaled to the square of the excitonic Bohr radius, Ob- 



III. TWO-DIMENSIONAL EXCITON WAVEFUNCTION 



In this section we calculate the wavefunction of an electron-hole bound state in the statically screened potential, Eq. 
@. We suppose that the interaction between the two particles depends only on the relative distance p = \r e — r^\, 
and we split the problem into two parts: the study of the relative motion of the two particles and the study of the 
motion of the center of mass which does not depend on the interaction. The total Hamiltonian can then be written as: 
Htot = H cm + H xe i- To study the relative motion problem we apply the variable phase method of scattering theory, 
neglecting the intrinsic spin effects. 

The Schrodinger equation for the radial wavefunction of the relative motion has the following form: 



dp 2 



1 d 
pdp 



,2 I 



U(P) 



p 



Rm.AP) = °> 



(6) 



for a given value of m , the azimuthal quantum number and where k 2 = 2m I E/ti 2 and U(p) = 2m T V s (p) j h 2 . 

For the bound states, the energy E is negative and we introduce the imaginary wavenumber k = in. As the 
potential vanishes at large distances, the solution of the radial equation Eq. (|16[) can be approximated for large p by 
the solution of the free Bcsscl equation, which is a linear combination of the modified Bessel functions of the first and 
second kind. Then, the solution of the radial Schrodinger, Eq. ([To]) , can be written as follows [301 ] : 
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(Kp) cos r) m -\ — K m (Kp)smr] m , (7) 

where I m (np) and K m (np) are the modified Bessel functions of the first and second kinds respectively while the phase 
shift r] m characterises their admixture and A m is the wavefunction amplitude. 

To solve the problem for all p, not just p — > oo, the phase shift r\ m and the amplitude A m are both considered not 
as constants but as explicit functions of p and k in the 2D formulation of the variable phase method (33j . Following 
Ref. [30|, we insert Eq. ([7]) into Eq. (J6j) and find that the phase shift satisfies the following first order, non- linear 
differential equation of the Ricatti type: 

d 7T „, s / , s , s 2 " 



, Vm, K {p) = -77 P ?7(p) x I m (np) cos7? m K (p) + - lf m («;p)sinr? m!re (p) 

Equation (|8|) is called the phase equation and should be solved with the boundary condition: 7y m-K (0) = 0, thus 
ensuring that the radial function does not diverge at p = 0. For the bound states the diverging solution vanishes, 
thus implying the asymptotic condition: 

lim Tj m ,K(p) = {v- l/2)7r, (8) 

p — >oo 

where v enumerates the bound states for a given m. The number of non-zero nodes of the radial wave function is 
given by v — 1 . 

Similarly the amplitude A m ^ K (p) satisfies the following equation: 

2 

(np) suvq m , K (p) K m (Kp) COS Tj 77l Ki 

-T- An, K (p) = A -m, K {p) \ X — rj mtK {p), (9) 

(kp) cosr] m , K (p) + - K m ( K p) Sin 37m, re (p) IJ 

TT 

which is coupled to the phase equation, Eq. ([5]), whose solution is computed first. 

The total bound exciton wavefunction is hence given by the product of the relative wavefunction 



ip m , K (p,ip) = I I m (Kp) cos r] mtK (p) + - K m (Kp) sin r] mtK (p) J x A m:K (p)e tmv , (10) 

and the exciton center of mass motion wavefunction characterized by the plane wave <^k cm (-R) = exp(— ik cal ■ R)/y/A., 
where R — (m c r c + m^r^) / (m c + m^) is the exciton center of mass coordinate, r G and the carrier coordinates , 
A is the surface area of the 2D system and fc cm = \/2ME cm /h the exciton center of mass momentum. Calculation of 
the scattering states is detailed in Ref. [3(| ■ 



IV. PHOTON EMISSION INDUCED BY EXCITON-CARRIER SCATTERING 



A. The scattering matrix elements 



In this section we are concerned with the scattering of Is excitons with free carriers: the excitons transfer their 
momenta and energies to the free carriers and reach the photon line where they recombine and thus emit a photon. 
Two possible spin configurations should be considered: depending on its total spin state an exciton may or may not 
radiatively recombine because of selection rules. In the X — e~ scattering process, the triplet state is defined with 
parallel electron spins and the singlet state with antiparallel electron spins (the same obviously applies for hole spins 
in the case of X — h scattering). If, for instance, the exciton spin is +1, the exciton is said to be bright and if it is 
+2, the exciton is dark. Therefore, exchange between the free carrier and the bound carrier changes the nature of the 
exciton in the singlet configuration, while the total exciton spin is conserved in the triplet configuration. 

The carriers interact via the screened potentials defined in Eq. Q and the matrix element V^ t , c = e,h, can be 
written as a sum (difference) of a direct term and an exchange term V^ h in the triplet (singlet) configuration: 
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Kf t (fc cm , fc 2 ) = V d f (fc cm ) ± V c f ch (k cm , fca). (11) 
As shown in Appendix A, in the case of exciton-electron scattering: 



T/cX _ e 2 i? 2 , K (7hfccm) -flp. B (7cfcc m ) 

Mir — o J ~ ' V i/! J 



and 



Kxch( fe cm,fc 2 ) = / — ! dq + i? , K (fc 2 +7hfccm) / — ; ; dq 

^0^1 \J q + 9s J q + q s 

- R QAk2 ) I ^fo+T^-q) J , (13) 

7 9 + 9s J 

where j c = m c /M. Note that the minus sign before the third term in the equation above may change the overall 
sign of the exchange term when it is computed as a function of k cia . To obtain these two terms for the exciton- hole 
scattering, one should simply swap the effective masses m e and where appropriate and change the overall sign. 

In Eqs. (fT2")) and (|13p. the quantities Rq and i?g respectively denote the Fourier transforms of the m — radial 
wavefunction in Eq. (|10[) and its square. Since the free carriers wavefunctions are taken as plane waves, the direct 
term V^f, c=e,h, is independent of the scattering free carrier momentum fc 2 and only depends on the energy of the 
bound state characterized by the wavcnumber k, the effective masses of the carriers m e and m^, and the initial exciton 
momentum fc cm . We also find that V^(k cm ) = -^(tm). 

B. Emission rates 

In this section we give the expressions of the exciton-carrier scattering contributions to the photon emission rates, 
i? e x(Sfi) and i?hx(^)- The expression for i? eX (M7), derived from Fermi's golden rule, is taken from Ref. 



with 



- 1 W ( (1 Jg^ +4< ) *> - « - + »- -Ml. <*) 

where 14cat(fccm, ^2) are the exciton-carrier scattering matrix elements given by Eqs. (|12[) and (|13|) . fe cm the exciton 
center of mass wavevector, and fc 2 the free electron wavevector. The exciton and free carrier distributions, Nk cia and 
N k2 , are: 

^ = ^-(l-«)nexp^^j, (16) 

and 

N k2 = ^— an exp -/? tH" ■ ( 17 ) 



The coefficient £ in Eq. (fT5|) ensures that there is no divergence when the photon energy Ml is equal to the gap energy 
E g . The exciton binding energy is E*, so the total energy of an exciton is: Ex = E g — E^ + h 2 k 2 m /2M . 



7 



TABLE I: Values of plasma screening parameter q s scaled to the excitonic Bohr radius, and the corresponding densities n in 
10 11 cm -2 and degrees of ionization a, for ZnSe and GaAs QWs at T — 300 K. The exciton binding energies, E£ , are given in 
excitonic Rydbergs. 



g s aB 


n 


GaAs 


a(n) 


/; 


ZnSe 


a(n) 


El 


0.10 


0.14 




0.93 


0.66 




0.41 


3.30 


0.32 


0.48 




0.86 


2.51 




0.34 


2.51 


1.0 


1.62 




0.78 


6.78 




0.39 


1.35 



To proceed with the calculations we combine Eqs. (|14lll7j) . which yields: 



Rcxifrty = -r~, n % , 9n9 ° > (1 - o)n exp \-p — (18) 

for the exciton-electron scattering. The function R^x(Ml) is formally identical to i? e x(^) in Eq. (fT5|) . The effective 
masses m e and mh have only to be changed where appropriate. The rates R c x{t>£l) depend explicitly on the degree 
of ionization, a, via the product a(l — a). This product is maximum for a = 1/2, which means that for given 
thermodynamic conditions on temperature and density, i? c x(^), as defined above, is greater when the plasma is 
composed of a mixture of excitons and free carriers in equal proportions. As mentioned in section II. B, the nature 
of the electron- hole plasma in ZnSe QWs is closer to a plasma with a = 1/2 than it is in GaAs QWs, and one 
should expect a more important contribution of exciton-carrier scattering to PL in ZnSe QWs. Note that very high 
proportion of excitons do not guarantee a great contribution to i? c x(^) either, as it would simply mean that there 
would not be enough free carriers to scatter with. The key point here is to have a close to 1/2. 

For numerical purposes, the discrete sums are approximated by 2D integrals: — > A/An 2 J dk. To evaluate the 
effects of exchange on the emission rates induced by exciton-carrier scattering, we first calculate R c x(fiQ), considering 
only the direct term of the scattering matrix elements. Note that, as shown in Appendix B, the numerical calculation 
of R c x(MY) is greatly simplified in this case since does not depend on (see Eq. (|B8[) L Inclusion of the exchange 
terms in the calculation of i? c x(^) requires the full numerical computation of Eq. (|18p . 



V. NUMERICAL RESULTS 



The materials parameters for GaAs-based quantum-wells are: m c = 0.067rno: m h = 0.46too ; where mo is the free 
electron mass, e r = 12, E g — 1.45 eV. For ZnSe we use: m c = 0.15mo, = 0.60mo, e r = 8.8, E g — 2.7 eV. As 
in Ref. 0] we take £ as 10 -3 . The ionization degree is evaluated at T = 300 K for various plasma densities given 
in Table 1. To study V^ t as a function of the exciton center of mass wavevector, fc cm , we fix the modulus of the 
free carriers wavevectors, fc2, evaluated from the mean velocity of the free carriers within the 300 K carrier gases. In 
GaAs QW, &2 = 2.32 a^ 1 and ki = 6.07 a^ 1 for the electron and the hole respectively. In ZnSe QW, these values are 
k 2 = 1.24 a~r} and k 2 = 2.48 On 1 . 



A. The scattering matrix elements 

To calculate the direct and exchange terms of the scattering matrix elements, Eqs. (ji"2"]) and (fT3")) . we need to 
compute first the relative motion part of the exciton wavefunction in the screened Coulomb potential from the 
numerical solution of Eqs. ([8]) and ([9]) obtained with the variable phase approach. In Fig. [2] the behavior of the direct 
terms of the scattering matrix elements versus A: cm is shown for two values of the plasma screening: q s ciB — 0.1 and 
q s aB — 1, for both materials. In all cases the qualitative behavior of is the same: V^f = for k cm — since the 
exciton cannot scatter from fc cm = to fc cm = 0; then, for small wavenumbers, which is negative for all values 
of A; cm , decreases and reaches a minimum for fc cm w Sa^, 1 , which means that the location of the minimum mostly 
depends on the ratio of the electron and hole effective masses, and less on the screening parameter in the low density 
regime. For exciton center of mass momenta larger than 3 a^ 1 , V^."^ e is a monotonically increasing function of fc cm 
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whose amplitude thus diminishes: the transfer of increasing large exciton center of mass momentum to its scattering 
partner is less likely. The matrix element is, as expected, greater for low screening, but not in a dramatic way: 
in the low density regime the screening parameter remains small enough not to have a siginificant impact on the 
amplitude of V^. The magnitude of is also higher in ZnSe than it is in GaAs, reflecting the greater strength 
of Coulomb interaction in ZnSe. Note that if the electron and hole effective masses were equal, the electron-electron 
and electron- hole contributions to would cancel exactly and the scattering amplitude would be zero for all fc cm - 
For a given plasma screening, the same comments apply to except that because of opposite signs, it is positive 

for all fc cm . The functions and VjSF reach their extrema for the same values of fc cm . Note that the extrema of 
V£fr, c = e,h, in both ZnSe-based and GaAs-based QWs, are reached for values of k cm that are very close: here the 
material dependent effective masses have a limited influence on the behavior of the function V^(k cm ). 



ZnSe QW 



GaAs QW 
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FIG. 2: Direct terms of the scattering matrix elements as functions of k C m, describing the Is exciton-carrier direct scattering 
within a screened Coulomb potential in ZnSe and GaAs quantum wells, in a low density regime at T — 300 K. For a given 
plasma screening is the exact opposite of V^. 



We turn now to the study of the full scattering matrix elements, V^at; including the exchange terms. For discussion, 
the value of is taken as the thermally averaged wavevector as introduced at the beginning of this section. A 
comparison of amplitudes of and V** h is depicted in Fig. [3J For a given value of k%, V^ h is much greater than 
V^. Exchange effects are therefore of importance even if the plasma is at room temperature. Another interesting 
observation is that unlike the direct term, the exchange term changes sign, thus reflecting the complex interplay of the 
Coulomb interactions between the three scattering partners. These results are consistent with the findings of Ramon 
et al [38j in the case of electron-exciton scattering. Note that the calculations of Ramon et al [38| were performed in 
the low temperature regime thus resulting in a greater amplitude of the exchange term at fc cm ~ 0. 

The total scattering matrix elements can thus be safely approximated by the exchange terms as shown in Fig.[U For 
a given plasma screening, we find that the scattering matrix elements, V^ t and V^ t , are not the exact opposite of each 
other, unlike the direct terms. This originates from the non-trivial dependence of V^ h on the free carriers wavevectors 
&2 and the effective masses as can be seen in Eq. (|13[) . For small and moderate exciton transferred momenta fc cm , 
the amplitudes of the scattering matrix elements V££ t are smaller than those of V^at simply because the exchange 
effects are less important with larger effective mass. However, as the value of fc cm increases the amplitudes become 
comparable and we observe that they decrease in a slower fashion for Vj££ t than for V^ t . Because of their smaller 
effective mass, the electrons acquire a larger kinetic energy for a given value of the transferred momentum than the 
holes, which as a consequence diminishes the exchange effects for the lighter quasiparticles. The interplay between 
effective mass and transferred momentum thus has a non-trivial influence on the exchange term of the scattering 
matrix elements. 



B. Contribution of exciton-carrier scattering to the emission rates 

The numerical evaluation of Eq. (|B8[) gives the Is exciton-carrier direct scattering contributions to the emission 
rates, shown in Fig.0 For a given material the behavior of R c x(tiQ) reflects its complex density dependence via the 
factor cv(n)(l — a(n))n 2 . Note that since the scattering matrix element does not depend dramatically on n it has a 
rather small influence on the emission rates: i? c x(^) is proportional to the square of and while the amplitudes 
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FIG. 4: Scattering matrix elements (approximated by the exchange terms) as functions of fc cm , describing the Is exciton-carrier 
scattering within a screened Coulomb potential in ZnSe and GaAs quantum wells, in a low density regime at T — 300 K. 



of decrease as the plasma density increases, the function i? c x(M7) is increasing significantly. It is also of interest 
to note that the magnitude of R^xiMY), induced by the direct term of the scattering matrix elements, is greater than 
R c x(fiD,) for comparable densities. This can be explained as follows: first, since the direct terms of the scattering 
matrix elements and have exactly the same magnitudes for a given carrier density, the X — e~ direct 
scattering channel is not favored despite the electron lighter effective mass; however, as can be seen from Eq. (|B8I) . 
R c x(Ml) increases with increased carrier effective mass m c (which can be found in the prefactor before the integral). 
Second, we checked that the function of photon energy, H£l, resulting from the computation of the integral in Eq. (|B8[) 
(without the prefactors) is greater in the case of X — h direct scattering near its maximum. This mass effect can be 
interpreted in terms of the two-dimensional density of states that is indeed greater for the hole gas than it is for the 
electron gas. Note that the emission spectra exhibit a thermal tail on the low energy side of the spectrum mirroring 
the free carriers distribution functions. 

To calculate the emission rates, including exchange effects in the scattering process, from Eq. (|18p . we should, in 
principle, take the carrier spins into account and hence study both the singlet and triplet configurations. However, as 
clearly shown in Fig. |4j the direct term is negligible compared to the exchange term. Therefore, the total scattering 
matrix elements can be very well approximated by the exchange term only in Eq. and the spin configuration 

becomes irrelevant here, as observed by Ramon et al [38j . In Fig. [51 we find that the inclusion of exchange effects 
in our model reverses the trend that we observed above: the magnitude of Rex(MY) has increased dramatically and 
is now greater than that of i?hx(^)- Exchange is therefore the mechanism that makes the X — e~ scattering a 
more efficient process than X — h scattering to lead an exciton to the radiative cone where photon emission occurs. 
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FIG. 5: Contribution of the direct Is exciton-carrier scattering to the emission rates evaluated for ZnSe and GaAs QWs in 
the nondegenerate regime at T = 300 K. The corresponding values of the densities for the screening parameters are given in 
Table [I] The rates of photon emission induced by the direct exciton-hole scattering dominate those induced by the direct 
exciton-electron scattering in both materials, by one order of magnitude. Note the differing vertical scales. 
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FIG. 6: Contribution of the Is exciton-carrier scattering to the emission rates evaluated for ZnSe and GaAs QWs in the 
nondegenerate regime at T — 300 K. The corresponding values of the densities for the screening parameters are given in 
Table [T] Exchange has reversed the previously observed trend: it is now exciton-electron scattering that yields the highest 
photon emission rates in both materials. Note the differing vertical scales. 



Although of importance, exchange effects are indeed reduced for the holes, which are heavier (4 times and 7 times the 
electron effective mass in ZnSe and GaAs respectively). Comparing results obtained for ZnSe and GaAs, it appears 
clearly that, for a given plasma screening, the contribution of exciton-carrier scattering to the emission rates at room 
temperature is more important in ZnSe than it is in GaAs, whatever the scattering channel is, X — e~ or X — h. 
This is first due to the fact that for a given plasma screening, the plasma density is higher in ZnSe QW than it is 
in GaAs QW; second, because of a smaller dielectric constant that yields a strong exciton binding energy in ZnSe 
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QWs, the exciton population remain important at low and moderate carrier densities (more than 50 %) whereas the 
fraction of excitons in GaAs QWs remains modest (less than 20 %) and hence yield a lower carrier-assisted excitonic 
contribution to photon emission. The nature of the mixed exciton/electron-hole plasma is therefore of importance. 

VI. SUMMARY AND CONCLUSION 

We have computed the exciton-electron and exciton-hole scattering matrix elements and discussed their properties. 
We have considered the fermionic nature of the carriers through exchange effects, which have been found to be 
very important even in a dilute mixed exciton/electron-hole plasma at room temperature. In these thermodynamic 
conditions the carrier-assisted radiative exciton recombination process has been investigated for GaAs-based and 
ZnSe-based quantum wells. For a given plasma screening, the emission rates due to elastic exciton-carrier scattering 
in ZnSe-based QWs are found to be one to two orders of magnitude greater than those in GaAs QWs. This is partly 
due to the fact that for the same value of plasma screening, the plasma density is greater in ZnSe QWs than it is in 
GaAs QWs, but more importantly, because of the lower ionization degree of the ZnSe QWs plasma. We have also 
found that in both systems the X — h scattering channel is a more efficient one in the case of direct scattering when 
the densities of electrons and holes are comparable. We attribute this essentially to the two-dimensional density of 
states which is greater for hole than for electron systems. The direct terms of the scattering matrix elements have 
been found to be negligible compared to the exchange terms. As a consequence, fermion exchange favors the X — e~ 
scattering channel for the carrier-assisted radiative recombination. 
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APPENDIX A: CALCULATION OF THE SCATTERING MATRIX ELEMENTS 

In this appendix we show how the two terms of the scattering matrix elements whose expression is given below, are 
calculated in the case of exciton-electron scattering. We write V^ t (fe cm , k 2 ) as follows: 

V°* t (k cm ,k 2 ) = \J dndradn, (^L +feom (r2)*S ire)ff (0,p lh )-4 +fecm (r 1 )*J )K5 (0 ) p 2h )) (Al) 

X (V s (pih) + V s (p 2 h) + V s (pi 2 )) (*O, K ,fc cm (i?l,Plh)0fe 2 (T , 2 ) - *0,«,fe om (-R2,P2h)0ki(ri)) , 

where R\ and R 2 are the exciton center of mass coordinates, r\ and r 2 the electrons coordinates and the hole 
coordinate. The relative distances are defined as follows: pih = |rjj — r±\, p 2 ^ = | T*h — ^"2 1 and p\ 2 = \ri — r 2 \- Note that 
the attractive electron-hole interaction potentials Ks(pih) and Kj(p 2 h) are negative and the repulsive electron-electron 
interaction potential V s (pi 2 ) is positive. Since we consider Is exciton states only, the value of the projection of the 
angular momentum, to, is zero and we remove the explicit ip dependence from the expression of ^o,K,k cm - The initial 
state of the free electron is the plane wave <f>k characterized by the wavevector k (= k\ or k 2 ), and its final state 
the plane wave 4>k+k cm characterized by k + fc cm , where fc cm is the exciton center of mass momentum that has been 
transfered during the scattering process. Equation (|A1|) is similar to the expression in Ref. [32| . 

Starting from Eq. (|A1[) . the first step consists of identifying the relevant physical terms: the direct and exchange 
terms, V£j£ and respectively. Calculations yield the following expressions: 

= J 4 2 +fe c >2)*J K i5 (0,pi h ) (V s ( P2h ) + V s ( Pl2 )) *o,«,fe om (i*i,Pih)^ 2 (r-2) dndr 2 dr h (A2) 

and 

V c f ch = - f 4 2+fccm (r 2 )*J iK 5 (0, Plh ) (V s ( Plh ) + V s ( P2h ) + V s ( Pl2 )) *o, K ,fe cm (i* 2 ,P2h)<Mn) dndradrh. (A3) 
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An additional term V^dd appears: 

Kdd = J <t>i 2+kc Jr2)^l tKji {0,pi h ) V s { Plh ) *o, K ,fc cm (^i,Pih)0fc 2 (r 2 ) drxdr 2 dr h , (A4) 

and before proceeding with the calculations of V?£ and V^ h , we show that I4dd = 0. Recalling the explicit expression 
of the exciton wavefunction, ^ m K & , in Section III, the additional term may be rewritten as: 



T4dd ^ 



We define the Fourier tranforms of the square of the relative motion part of the excitonic wavefunction, Eq. (fT7|) : 

= Y.Rl,Mi) eW*-^ = ^J RIM e^^-^dqx, (A6) 
and of the scattering potential V s (p-)- 

V s ( Plh ) = 5^.(9) e^-^ = Ay V s (q) e^-^dq, (A7) 



which we include into Eq. (|A5|) to find: 



KkM = J^JJ e^^e- i ^*™ + « 1 +«> ri e- i ^ fc »-« 1 -«^J^ iB (g 1 ) V,(qr) dndradrhcUfcdqr, (A8) 

where 7 G = m c /M and 7h = rrv^/M. Now, considering the identity J e I<3 ' r dr = 47r 2 <5(q), where 5 denotes the Dirac 
function, the above integral reads: 



Kdd = J J 5hck cm + qi+q) S(-~/ h k cm + q)e lk ^ r2 Rl K {q)V s {q) dr 2 dqidq, (A9) 
and, since J 5 (a — x)S(b — x)f(x) dx — when a ^ b, we finally find that V a dd = 0. 

1. The exchange term 

We turn now to the exchange term V^. h that we artificially decompose into three terms: V^. h = V^. h 1 + V^ h 2 + 

T/oX . ' '' ' 

cxch,3 ' 

Kxch.i = - / 4 2 +fe cm (^)*J, Ki o^"° lh ) ^ lh ) *o,K, fccm (i?2,P2h)0 fcl (r 1 ) dndr 2 dr hl (A10) 
K^h, 2 = ~y 0l 2+fccm (r- 2 )vI/J KO -(O,p lu ) y s ( P2U ) *o, K , fccm (i? 2 ,p 2h )0 fcl (n) dndradrt, (All) 



and 



K^h,3 = - J 4 2+fccm (n)^ K O -(0,p lh ) K(pi 2 ) *o, K , fccm (fi 2 ,p 2 h)0 fcl (n) dndradrh. (A12) 



In the same fashion as above we define the Fourier tranforms: 
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^(P2h) = £ F-(q) e Mrh-r») = _^ J V-(q) e^-^dq, (A15) 



and 



^s(pi2)^K + (g) e M-*-rO = A. [ V + {q) e ^(r 2 - ri ) 



(A16) 



where the + and — signs denote the repulsive and attractive potentials. Inserting the relevant Fourier transforms of 
Eqs. (KE\i , (fATf . (|A13tiA16|) into Eqs. (|A10|) , (jAlip and (|A12|) and performing the same type of calculations as for 



V s f a t we find: 



and 



= / - q)RoA k 2 + 7hfccm + q)V s {q) dq, (A17) 

Kx*h,2 = -ARo,«(fe 2 + Thfecm) /" - q)V s -{q) dq, (A18) 

KSh,3 = -^Ro,«(fe) / ^o, K (fc 2 + 7hfccm - q)K+(g) dq. (A19) 

The exchange term, V^. h , depends on the properties of each scattering partner: exciton binding energy, kinetic 
energies, effective masses of the electrons, hole and exciton. To obtain the exchange term for the exciton-hole 
scattering, one simply needs to swap the effective masses m c and where appropriate and change the overall sign. 

2. Direct term 

Finally, we show the main steps of the calculation of the direct term starting from Eq. (|A2[) . We write — 
V^ T + V^7 r and we need only to explicit calculations for V^ T as V d t has a similar structure. V^ r is the first of the two 
terms in Eq. (fA"2|) : 

= J 0i a+)W (ra)*5 )Biff (ff J Pih)V B (/o^)*o J »,h < »(-Ri J Pih)0fc a (^) dndr 2 dr h (A20) 

Using the explicit expressions of the electron plane wave and exciton wavefunction of Section III as well as the 
definitions of the Fourier tranforms of Eqs. (|A6[) and (|A15D . we find: 

U d T r = — ^ / e «( fc ™+9)-^ e i(n-i,h cm )- ri e -i( 9 +qi+7hfe om )-r h ^2 dridr 2 dr h dqdqi. (A21) 

167T J 

This 10-dimensional integral can be simplified using again the definition the 5-function given above: 
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Vdh = ^,«(7efccm)V-(fe cm ). (A22) 

Doing the same calculations with V^ at , Eq. (|A2|) can be reduced to the following expression in the Fourier space: 

y d f (fe cm , K )=A (^, B (7efccm)V7(fc C m) + R$, K (-^cm) K + ■ (A23) 

Since the Fourier transforms of the interaction potential and the squared wavefunctions are radial functions, V^(fc cm ) 
only depends on the modulus of the vector fc cm . Combining Eq. (|A23[) with the Fourier transform of Eq. (|2"|). we 
obtain the final expression of the direct term of the exciton— electron scattering matrix elements, which reads 



V^(k cm , k) = - ^ °' kA '-. (A24) 

Note that unlike the exchange term, the direct term only depends on the exciton properties: kinetic energy, binding 
energy and effective masses. We find that V^{k cra ) is the exact opposite of V^(k cm ). This is due to the fact that 
both free carriers wavefunctions are assumed to be plane waves. 



APPENDIX B: CONTRIBUTION OF THE EXCITON-ELECTRON DIRECT SCATTERING TO 

PHOTON EMISSION 



To calculate the contribution of the exciton-electron direct scattering to photon emission, i? c x, we start from the 
evaluation of the discrete sum over all the vectors fc 2 in Eq. (fl"8]) . In this case V^ t (k cnl , fc 2 ) reduces to V^(k cm ) so 
that the scattering matrix element can be taken out of the sum over the vectors fo. To proceed, one can approximate 
the sum S: 



5>» ("' 3 3) 5 ( E » - E * + w - hn - £«■ + 2fc » ' fe > I • < B1 > 



by a 2D integral: 



The argument of the S function in Eq. (|B2[) . is of the form X — Y cos9, where 9 is the angle between the vectors 
fc cm and k2- The integral over 9 , I[5], can hence be calculated analytically: 



/■27T 

L ' 



5{X - Fcos 



if \X\ < \Y\. We find: 



|ysin[cos _1 X/y]| 



(B3) 



I [8] 



^fc cm fc 2 sin 



h 2 k cm k 2 /m e 



(B4) 



where X = Ex — Ml — h 2 k 2 m /2m e , Y = f?k cm kijm ei and with the above condition on X and Y, the limits of the 



integral Eq. JBSJ) are: kf m = [E x - HQ ■ 
Considering the identity |sin[cos _1 0] | 
to: 



h 2 k 2 nl /2m c ) m c /h 2 k cm and fc™ ax 



yfl - cos 2 [cos- 1 6] = Vl - B 2 and combining Eqs. (fB2|) and (|B4| lead 
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27T 2 Jfemin j n 4 k 2 .2 / t2.2 \ 2 



(B5) 



2to p 



Note that the analytical calculation would stop at this stage if the full scattering matrix element V r sca t(A: cm , k 2 ) were 
considered, since it would explicitly appear in the above integral. In the case of direct scattering, Kicat only depends 
on k cm and hence we may proceed with the analytical calculations as follows. 

With two successive changes of variables: K% — fcf an d K' 2 — i^/fc™" , Eq. (jB5|) becomes: 

Am e kf n f°° \ 2m 



S = 9 % 2 / i — '- AK' 2 . (B6) 



Considering the identity: e ^ x (x — 1) 1 I' 2 Ax = \/ttJJx e M , we find: 

4ir 2 H 2 k cm \ /3h 2 J XP \ ' 2m c J ' ( ' 

Finally, inserting Eq. (|B7|1 into Eq. (fT8|) and approximating the discrete sum over all the vectors k cm by an integral 
lead to the expression of the emission rate, i? e x(M7): 



k(1 - a)iV 2 V27rm e /3 3 4n£Ml/E s 



exp — P 



2 1.2 



h z k 



M 



+ 



(1 - h 2 VL 2 /E 2 y +4tt£ 



(B8) 



2M 2rfc; 



2 ,,2 
cm 



E v -E(: -Ml- 



rrih h 2 k 2 n 
™ 2M 



\V s ll(k cm )\ 2 dk c 



The expression of Rhx(MY) is simply obtained by swapping the effective masses to g and m-^ where appropriate. 
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